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Abstract. Optical microcavities are open billiards for light in which electromag- 
netic waves can, however, be confined by total internal reflection at dielectric bound- 
aries. These resonators enrich the class of model systems in the field of quantum 
chaos and are an ideal testing ground for the correspondence of ray and wave dy- 
namics that, typically, is taken for granted. Using phase-space methods we show 
that this assumption has to be corrected towards the long-wavelength limit. Gener- 
alizing the concept of Husimi functions to dielectric interfaces, we find that curved 
interfaces require a semiclassical correction of Fresnel's law due to an interference 
effect called Goos-Hanchen shift. It is accompanied by the so-called Fresnel filtering 
which, in turn, corrects Snell's law. These two contributions are especially impor- 
tant near the critical angle. They are of similar magnitude and correspond to ray 
displacements in independent phase-space directions that can be incorporated in 
an adjusted reflection law. We show that deviations from ray-wave correspondence 
can be straightforwardly understood with the resulting adjusted reflection law and 
discuss its consequences for the phase-space dynamics in optical billiards. 



1 Introduction 

Mesoscopic and nanoscopic systems [1 , ranging from quantum dots and 
nanoparticles to carbon nanotubes and most recently to graphene [2], still 
receive growing attention. These systems are characterized by a phase- 
coherence length that is larger than the system size such that interference 
effects play a crucial role and a quantum mechanical description is in order. 
Often, however, a semiclassical description is sufficient in order to explain 
the observations: mesoscopic systems, with typical sizes on the micrometer 
scale, are truly in the middle (greek messi) between the microscopic and the 
macroscopic world. The motivation for their investigation comes from both 
the application-oriented as well as from the basic-research side. Examples are 
the quest for the miniaturisation of electronic and optical devices or the chal- 
lenge to observe many-body effects known from bulk metals in finite systems, 
a prime example here is the observation of the Kondo effect in quantum dots 
[3] . Another inspiration comes from the field of quantum chaos [4] where the 
focus lies on the dependence of physical observables on the system geometry. 
More precisely, the crucial property is whether the dynamics of the under- 
lying classical system is chaotic or integrable, corresponding to a chaotic or 
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integrable (or, in the most generic cases, a mixed) phase space. The first 
observation of such a sensitivity was in the magneto conductance through 
circular (integrable) and stadium-shaped (chaotic) structures that shows a 
triangular and Lorentzian coherent-backscattering signature, respectively, in 
the magnetoconductance [5]. 

The objective of the present paper is to give an overview over some of the 
recent work on quantum chaos and semiclassical aspects in optical mesoscopic 
systems which we introduce in the following Section. We will then address 
a number of deviations from ray-wave correspondence and show how they 
can be explained when correcting the usual specular reflection law for light 
by semiclassical effects important on the mesoscopic scale. We end the paper 
with a discussion of the implications of such an adjusted reflection law. 

1.1 Quantum dots and optical microcavities: Billiards for 
electrons and light 

Quantum dots realized in semiconductor heterostructures are often the first 
systems that come into mind in the mesoscopic context. Electrons are con- 
veniently trapped and manipulated by various gate electrodes. Especially 
many-electron ballistic quantum dots are often described as billiards with 
hard walls andused as model systems in the field of quantum chaos [5]. 

Another class of experimental and theoretical quantum-chaotic model sys- 
tems are optical microcavities - billiards for light instead of electrons. This 
is possible because of the analogy between the Schrodinger and Hclmholtz 
equations for electrons and light, respectively, that holds in two dimensions 
and for the so-called TM (transverse magnetic) polarization direction (where 
the magnetic field lies in the resonator plane). One fundamental difference is, 
however, the confinement mechanism: For light, there is no charge to manip- 
ulate, and gate-voltage based confinement has to be replaced. The concept 
used instead is that of total internal reflection at optical interfaces with dif- 
ferent refractive indices n\ and 112- Total internal reflection occurs at the 
optically thinner medium (e.g., when going from glass with refractive index 
m = n « 1.5 to air with ri2 = 1) for angles of incidence x (measured to the 
boundary normal) larger than the critical angle \ c — arcsin(l/n). The limit 
n — > oo corresponds to that of a closed system with hard walls. This implies 
in turn that generic optical systems are open systems where the openness is 
related to, and defined by, the possibility of refractive escape. In other words: 
Optical microcavities are ideally suited to theoretically as well as experimen- 
tally study quantum (wave) chaos in open systems. Note that the openness 
is not induced by leads and that we are not interested in transport through 
the cavity. Rather, the openness exists everywhere along the boundary which 
is reflected in mixed boundary conditions for the Helmholtz equation with 
nonvanishing wavefunction and nonvanishing derivative. 




Fig. 1. Ray and wave description of the dielectric disk, (a) Whispering-gallery type 
ray trajectory, (b) Resonances of the closed (left half, corresponding to a wavenum- 
ber kR = 9.761) and dielectric (right half, n = 1.54, nkR = 11.428 - 0.254i) disk. 
The openness of optical systems considerably changes the wave intensity outside 
the cavity (evanescent wave, or refractive escape), (c) Fresnel's law with devia- 
tions from the naive ray- wave correspondence clearly visible in the region of critical 
incidence. See text for details. 



1.2 Ray- wave correspondence!? 

One paradigm intimately related to the field of quantum chaos is the 
quantum-classical, in the case of optical microcavities, the wave-ray corre- 
spondence. It is usually taken for granted and one of the footings of our 
present understanding of the relation between the classical and the quantum 
world. We shall see below that, although useful and extremely simple to im- 
plement numerically, the ray picture is not able to comprehensively explain 
all the results found for optical microresonators. A detailed list of observed 
deviations will be given in the following section. Note that these deviations 
are related to the fact that the wavelength A of the electromagnetic field is 
smaller, but still comparable to the cavity size R (and, therefore, the true ray 
limit A — > is not fulfilled): In typical experimental setups based on semicon- 
ductor heterostructures (n = 3.3), the cavity size is of the order R ~ 50^m 
and infrared light with A ~ 850 nm is used. In numerical wave simulations, 
values nkR exceeding, say, 250, are very hard to reach at present. 

1.3 Description of optical microcavities 

At this point, a few words concerning the description of optical microcavities 
in the ray and wave picture are in order. 

Ray picture: The ray optics description of optical microcavities is based on 
ray tracing simulations. The trajectory is determined by assuming specular 
reflections at the resonator walls. In addition, one introduces a variable that 
monitors the (decaying) intensity of the light ray given by Fresnel's law. It is 
even possible to account for the interference of rays [6] . Predictions made for 
the experimentally accessible far-field intensity are now mostly based on the 
refractive escape of rays taken from the steady-probability distribution [7] . 
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ThePoincare surface of section summarizes the information about the reflec- 
tion points at the (outer) resonator boundary in reduced phase space given 
by a spatial variable parametrizing the position along the boundary (such as 
the polar angle <f> or the arclength s) and the angular momentum sinx, see 
Figs. Ufa) and[2jc). The condition for total internal reflection, | sinxl > l/n, 
is violated in the so-called leaky (forbidden) region —1/n < sin^ < 1/ti 
[marked by dashed-dotted lines in Fig.[2jc)]. Trajectories hitting this phase- 
space region will (more or less) easily escape the cavity by refraction. Periodic 
orbits with stable islands in the leaky region are considered to be not popu- 
lated by cavity modes. 

Wave Picture: The objective is to compute the resonances (or quasibound 
states) of the cavity. The by far most popular approach and a versatile tool 
is the boundary element method [8] that gives the resonances directly in the 
complex plane. The imaginary part of the dimensionless complex wavenumber 
kR (with k = 2tt/X and R the radius of curvature) in free space defines the life 
time and the Q-factor of the resonance, Q — l/[2 Im(fc-R)]. An example of a 
resonant wave pattern is shown in FigQJb) where that of a closed system (left 
half) is compared with that of an open, optical system with refractive index 
n = 1.54 (right half). Another approach is to use an 5-matrix method [9] ITD] 
that describes the resonator as an open system being "probed" from outside 
with (test) plane waves. The position of resonances can be read-off from the 
Wigner delay time. This approach is straightforwardly implemcntedfor (at 
least partially) rotationally invariant systems such as the annular billiard, 
see Fig.CJb). 

Husimi functions at dielectric interfaces: The mapping of the resonance 
wave pattern to phase space is realized by means of the Husimi function [11] : 
Simply speaking, this function measures the overlap of the resonance wave 
function with a minimal-uncertainty wavepacket centered around a certain 
position 0OjXo i n phase space. As illustrated for the example of the annular 
billiard in Fig. [3J Husimi functions are a particularly useful tool to study 
ray-wave correspondence. We point out that in order to do so, the concept 
of Husimi functions has first to be generalized to dielectric interfaces [TTj . 
In hard-wall systems with Dirichlet boundary conditionsthe wave function 
vanishes at the system's boundary, and the Husimi function can be defined 
based on the wave function derivative alone. At optical interfaces, however, 
both the wave function and its derivative are non-zero and it is not clear at 
all which of the two should be used to define the Husimi function. Moreover, 
there are now four rays (incident and outgoing on either side of the interface) , 
and the existence of four corresponding Husimi functions would certainly be 
desirable and support the concept of ray-wave correspondence. As a matter 
of fact, both issues can (solely) be solved simultaneously as shown in Ref. [IT] 
to whichwe refer the reader for details. 

Active (lasing) microcavities: From an application-oriented point of view, 
the most interesting application of optical microcavities is that as microlasers. 
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Fig. 2. Ray-wave correspondence in the dielectric annular billiard (with refractive 
indices n a = 3 and m = 6 of the annular region and the inner disk), (a) Stable 
periodic orbit typical for the geometry shown, (b) (One) corresponding resonance 
wave pattern, (c) Phase-space representation in terms of the Poincare surface of 
section (green/light crosses) and the (inside incident) Husimi function (maxima 
in red/dark regions) that reveal deviations from ray- wave correspondence on the 
quantitative level. The dashed rectangle marks a region where violation of forward- 
backwardsymmetry (sinx — > ™ sinx) is especially evident. The leakyregion, char- 
acterized by refractive escape of rays, is located in between the dashed-dotted lines. 



This requires the presence of an active material that allows for lasing oper- 
ation. Then, in addition to the nonlinearities originating from the (chaotic) 
resonator geometry those from the lasing operation are important and de- 
termine the behaviour. Although a description of the active material can, to 
a certain degree, be achieved using a complex refractive index n |12j . appli- 
cation of the Schrodinger-Bloch model [13] is in order. We will not further 
consider the wave description of active cavities in the present paper. We 
point out that, interestingly, the ray picture (where no active medium can 
be accounted for) may provide a very reasonable description of the far-field 
characteristics even for lasing microcavities [14]. Given the goal of building 
microlasers with directional emission (which comes close to being the holy 
grail in this field at present), ray simulations have proven to be a valuable 
tool even away from the ray limit for A < R [15j . 



2 Deviations from ray-wave correspondence 

In the following, we list a number of deviations from ray-wave correspondence 
that have accumulated over the past, say 10, years. Though each individual 
fact might be considered just as a slight mismatch, in its entirety these ob- 
servations suggest that the ray description needs to be corrected away from 
the exact limit. The necessary corrections will be discussed in the following 
section, and we finish the paper with an outlook about further consequences 
of an adjusted ray model. 
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2.1 Delayed onset of total internal reflection in Presnel's law 

One nice example illustrating the interplay between the ray and the wave 
description is the dielectric disk. On the ray side, the rotational invariance 
of the system conserves the angle of incidence x- If X is larger than the 
critical angle, light is confined in a so-called whispering-gallery (WG) orbit, 
cf. Fig. [BJa). On the wave side, this is reflected in an azimuthal quantum 
number to (e.g., to = 3 for the WG mode in Fig. [Tfb); here, the radial 
quantum number p = 2) that further characterizes resonances with complex 
wavenumbers kR. There is a one-to-one relation between the ray and wave 
picture quantities (which are% and the Fresnel reflection coefficient Rp(x) on 
the ray, and m and kR on the wave side), namely [51 110) 
vn 

S{nX= nMkR)^ (1) 
Rp = exp(4nlm(fci?) cosx) ■ (2) 

This translation works well, except near the critical angle where the on- 
set of total internal reflection is significantly delayed in the wave picture, see 
Fig- [He). The green/light (solid) curve shows the Fresnel law (for TM polar- 
ization) with Rf = 1 for x > x c ~ arcsin0.65. The squares correspond to the 
Fresnel reflection coefficient, Eq. for resonances with Re(kR) w 50. Note 
that a closed analytical expression for Rp (dashed line) was derived and is 
given in Ref . [16] . 

The origin of the deviations near Xc must be, and is, the curvature of the 
dielectric interface. The question arises (especially when taking a quantum- 
chaos inspired point of view) whether a semiclassical explanation of this devi- 
ation is possible, i.e., one that is based on the ray picture but takes correction- 
soriginating in the wave character of light into account. That this is indeed 
possible can be seen by the purple/dark (solid) line for a ray picture com- 
pleted by such a wave (interference) correction known as the Goos-Hanchen 
shift pTl [18] that then closely follows the squares [16] . We will discuss this 
effect in detail in Section [3] 

2.2 Correspondence of orbits and resonances in configuration and 
phase space: Qualitatively, not quantitatively 

In Fig. [5] ray-wave correspondence is illustrated for the example of the di- 
electric annular billiard. For the geometry chosen, the trajectory/resonance 
shown in Fig. 0(a) in the ray, and in Fig.[^b) in the wave description belongs 
to one typical family of stable orbits (or resonances). The similarity between 
the two patterns is evident and certainly supports the concept of ray-wave 
correspondence. When going from configuration to phase space, however, the 
agreement becomes somewhat less convincing, cf. Fig. [2fc) . The larger red 
(dark) maxima of the (incoming) Husimi function coincide reasonably well 
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with the Poincare signature of the orbit (green crosses), but the smaller max- 
ima (corresponding to the right reflection points closer to the constriction) 
clearly deviate from the ray signature, which is at least partially related to 
the Goos-Hanchen shift mentioned above. Such a behaviour of qualitative, 
but not quantitative ray-wave correspondence is the typical case. We will 
investigate the reasons and mechanisms leading to these deviations below. 

2.3 Husimi functions reveal violation of 
forward-backwardsymmetry 

One intrinsic property of the ray picture is the principle of reversibility of ray 
trajectories - in other words, forward-backward or time-reversal symmetry. 
The corresponding symmetry operation, sinx — > — sin%, is strictly obeyed 
in the ray picture. This is, however, not the case in the wave description as 
can be seen in the Husimi function, Fig. ^c), especially in the area marked 
by the dashed rectangle. There are numerous examples of such a behaviour, 
e.g. [19] . that was recognized in a number of cases but could not be fully 
understood. Usually, one was content that time-reversal and spatial (mirror 
axis of the billiard, e.g., x-axis in the annular billiard with the symmetry 
operation <f> — > 2tt — cf>) symmetries together were obeyed. We shall see in the 
next section that another semiclassical effect, the so-called Fresnel filtering, 
is responsible for the violation of time- reversal symmetry and implies, when 
taken into account in an corrected ray picture, non-Hamiltonian dynamics in 
optical microcavities. 

2.4 Existence of regular modes in chaotic systems 

Optical microresonators studied recently in quite some detail both in theory 
[2Ti] and in experiment [5T] are cavities with spiral shape, r(<j>) = R(l+e<j>/2Tr), 
where e measures the radial offset. In terms of classical dynamics, the spiral 
billiard is (for all purposes of a physicist) chaotic. The more surprising was the 
observation of predominantly regular orbits of triangular and star-like shape 
reported in Ref. [20] . There are strictly no periodic ray orbits correspond- 
ing to these patterns, neither stable (those are missing in chaotic systems) 
nor unstable (that could become visible as so-called scarred resonances) . The 
observed regular orbits were named quasiscars, and vicinity of the angle of 
incidence to the critical angle was suspected to be crucial (indeed, triangu- 
lar orbits were found in cavities with n = 2, the star-like type for n = 3). 
Note that the existence of these regular orbits suggests that the system ap- 
parently possesses rotational symmetry. Clearly, this cannot be the case in 
a classical description of spiral billiard. The question arises what mechanism 
re-establishes the rotational invariance in the wave picture. We shall see be- 
low that it is again the Fresnel filtering effect that we will explain in detail 
in Sectional 
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Fig. 3. Reflection of (a) a light ray and (b) a Gaussian beam at a curved interface. 
Clearly, the ray picture prediction fails towards the wave regime A < R. The reason 
are Goos-Hanchen shift and Fresnel filtering, see also Fig. [5] that vanish only in the 
pure ray limit A/R — > 0. 

2.5 Violation of Presnel's and Snell's law at curved interfaces 

Eventually, we consider the (single) reflection of a Gaussian test ray at a 
curved interface. For convenience we choose as a cavity an air hole of radius 
R in a glass matrix, or equivalently, consider a hole with refractive index 
n = 0.66 < 1 in air. (Note that this implies concave instead of the convex 
curvature considered so far.) In Fig. [3] the ray and wave result are shown 
for the specific case of angle of incidence xo = 45° (controlled via the im- 
pact parameter s). Accordingly, the reflected light ray will leave the cavity 
vertically. This is, however, not the case in the wave description: Here, the 
reflected beam leaves the cavity under a larger angle, see Fig. [3^b). 

3 Correcting ray optics by wave effects: Goos-Hanchen 
shift and Fresnel filtering 

At this point a systematic study of the effects causing the deviations from 
ray-wave correspondence is in order. Inspired by the doubtless advantages of 
the ray model, such as its easy implementationand its conceptual success, our 
objective will be to identify semiclassical correctionssuch that the resulting 
(adjusted) ray model can quantitatively better capturethe wave properties 
of the system. Following the last example in the previous Section where de- 
viations between the ray and wave behaviour are clearly visible in a single, 
near-critical reflection, we analyze this situation in some more detail |22j . Our 
means of choice are Husimi functions at dielectric interfaces. Incident and 
outgoing Husimi functions for the reflection of a Gaussian beam at a cruved 
interface, cf. Fig. [3J are shown in Fig. [5] The intersection of the green/light 
(dashed) lines marks the position of the maxima as expected from the ray 
picture. 

In the near (and similarly in the far) field, the incident Husimi function 
exactly coincides with the ray picture expectation:This choice defines our 
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(a) Husimi near field (b) Husimi far field 




Fig. 4. Husimi functions for the reflection of a Gaussian light beam at an circular 
inclusion with smaller refractive index as shown in Fig. [3] The parameters of the 
beam correspond to a light ray incident under \o = 40°, the critical angle is Xc = 
41.75°. Shown are the Husimi functions of the incident and reflected beam in the 
near and in the far field, see text and [22] for details. The intersection points of the 
dashed green/light lines indicate the ray model expectations. Clearly, the signature 
of the outgoing Husimi function deviates in both phase-space directions from it due 
to Goos-Hanchen shift Z\0gh and Fresnel filtering A\f- 

initial conditions for the incident Gaussian beam (in Fig. 2] the angle of inci- 
dence is xo = 40°). Deviations become visible when looking at the outgoing 
Husimi function: The maximum deviates from the ray model prediction in 
two independent directions in phase space, marked by the arrows in the inset 
of Fig. [4] The shift Z\^gh in (^-direction is known as the Goos-Hanchen shift 
|17l I18j , and the shift A\f in x~direction has been termed Fresnel filtering 
[23] . Both effects are schematically illustrated in Fig. \5j[b) . Their magnitude 
depends on the wavenumber and is typically several degrees [22] , with the 
Goos-Hanchen correction being the larger. 

Crucial for the understanding of both effects is to realize that in optical 
microcavities each light ray is actually a light beam, i.e., composed of light 
rays with similar but not exactly equal angles of incidence. This becomes 
immediately clear when recalling that in the mesoscopic regime a light ray 
assumes a transversal extension of the order A. At a reflection point, thecor- 
responding angles of incidence will acquire a certain distributionbecause of 
the interface curvature [16j . In addition, an electromagnetic wave, e.g., a 
Gaussian beam, always contains a range of angles of incidence. 

Goos and Hanchen showed in a nice experiment in 1947 |17j that, for 
angles of incidence larger than the critical angle, this leads to an interference 
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Fig. 5. Goos-Hanchen shift (GH) and Fresnel filtering at (a) a planar and (b) a 
curved interface. At planar interfaces, a slight variation in the angles of incidence 
around xo gives rise to the GH shift zgh, a lateral shift of the reflected beam in 
the regime of total internal reflection. At curved interfaces, the GH shift Z\0gh 
is directly observable as a change in the far field emission direction. In addition, 
Fresnel filtering A\f leads to non-specular reflection and a deviation from Snell's 
law that further corrects the far-field emission characteristics. 

effect that results in the lateral shift (of the order A) of the beam along a 
planar interfacebefore it is reflected. This is illustrated in Fig. EJa). The re- 
flection can be thought of as to take place at an effective interface somewhat 
inside the optically thinner material. Whereas the angle of incidence is the 
same at the real and the effective interface for planar interfaces, this is not 
true at the curved optical boundary of microcavities [16] : The angle of inci- 
dence is smaller (larger) at convex (concave) boundaries. Applying Fresnel's 
law to this effective angle of incidence can quantitatively explain the devia- 
tions in the Fresnel reflections coefficient discussed above, cf. the purple line 
in Fig. [TJc) and Ref. [TB] for details. 

Fresnel filtering is even more classically to explain than the Goos-Hanchen 
shift as there is no underlying interference effect. In a collection of rays with 
angles of incidence around the critical angle, the rays with the largest angles 
will already be totally reflected whereas subcritical rays are still refracted. 
This can be seen in the signatures of the outgoing Husimi functions in Fig.|4j 
The faint signature below the horizontal dashed green/light line corresponds 
to the refracted beam (that leaves the cavity when hitting the boundary the 
next time). It is indeed composed of angles x < Xc ~ 41.76°. In turn, the 
Husimi signature of the reflected beam is shifted above thedashed line that 
defines the position in the ray model. (Note that the opposite shifts of the 
reflected and refracted beam ensure conservation of angular momentum.) In 
other words, the maxima of the incident (Gaussian) beam and the reflected 
(asymmetrically perturbed) beam do not coincide: Around critical incidence, 
the angle of the reflected beam is always largerby A\F,ci- also Fig.JHb), and 
the law of specular reflection, and consequently forward-backward symmetry, 
is violated. 
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Goos-Hanchen shift and Fresnel filtering have a common origin - the 
beam rather than a pure ray nature of electromagnetic waves - but act in 
orthogonal directions in phase space and, therefore, cannot be comprised 
in one and the same correction to the ray model. On the other hand, they 
also exhaust the number of possible corrections because there are no more 
independent directions in phase space. We shall see in the last Section that the 
different nature of Goos-Hanchen shift and Fresnel filtering manifests itself in 
strikingly different effects on the dynamics in optical billiards described with 
an adjusted ray model that takes the above-discussed non-specular reflection 
near critical incidence into account [24] . 

4 Outlook: Non-Hamiltonian dynamics in 
quantum-chaotic model systems 

Given the increasing activity in the field of optical microcavities and quantum 
chaos over the past years, the question arises whether effects of such a non- 
specular reflection law have not been observed before. We already mentioned 
that the deviations in Fresnel's law, cf. Fig. |Tfc), can be fully understood 
with the Goos-Hanchen shift. It can also qualitatively explain the differences 
between ray orbits and resonance patterns, Fig. [2[ via an adjustment at the 
reflection points with the steeper (near critical) angle of incidence. In fact, all 
deviations from ray- wave correspondence discussed in Sec.[2]can be addressed 
based on Goos-Hanchen shift (that is important for all angles of incidence 
X > Xc) and Fresnel filtering (that is important especially for x ~ Xcj- The 
two remaining examples, the existence of regular orbits in spiral microcavities 
and the breaking of time reversal symmetry in the Husimi functions, can be 
explained by Fresnel filtering [24] , 

In the spiral cavity, the filtering correction Axf in the outgoing beam 
re-establishes the conservation of angular momentum as it counteracts the 
change in curvature that decreases x- Therefore, regular orbits (similar to 
those in the disk) may exist again [24] . They are unstable and, consequently, 
host true scars. In general, a finite Axf destroys the principle of ray path 
reversibility. This is easiest seen when considering a ray with xo ~ Xc that 
leaves then under an angle xo + ^Xf > Xc- Tracing its trajectory in oppo- 
site direction will yield a (nearly) zero filtering correction (such a situation 
can easily be constructed), and the reflected ray does not coincide with the 
original ray. It is precisely this type of mechanism that causes the observed 
loss of time- reversal symmetry in the Husimi functions. 

Most remarkably, a billiard dynamics based on the adjusted reflection law 
leads to non-Hamiltonian dynamics [24] . Responsible is the filtering correc- 
tion that causes deviations of the Jacobian matrix determinant from unity, 
see [24] for details. The origin is the openness of the billiard, effectively de- 
scribed by the adjusted reflection law. Note that dissipative as well as attrac- 
tive dynamics with the formation of repellors and attractors, respectively, is 
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possible [24] . We are optimistic that more signatures of non-Hamiltonian ray 
dynamics will be identified soonand that quantum chaos in open systems will 
remain a fascinating research topic in the future. 
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